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We report density-functional-theory calculations of the electronic band structures and optical ab- 
sorption spectra of two-dimensional crystals of Ga2X2 (X=S, Se, and Te). Our calculations show 
that all three two-dimensional materials are dynamically stable indirect-band-gap semiconductors 
with a Mexican-hat dispersion of holes near the top of the valence band. We predict the exis- 
tence of Lifshitz transitions — changes in the Fermi-surface topology of hole-doped Ga2X2 — at hole 
concentrations ns = 7.96 x lO" cm"^, nse = 6.13 x lO" cm"^, and tite = 3.54 x lO" cm"^. 

PACS numbers: 73.61.Ga, 78.66.Hf, 73.20.-r, 71. 15. Mb 



I. INTRODUCTION 

Two-dimensional semiconductors have been studied in- 
tensively in the past few decades as researchers have 
tried to find materials to complement or even replace 
silicon in electronics. A recent trend in the hunt 
for new two-dimensional systems consists of the iso- 
lation and study of atomically thin sheets of layered 
materials. Honeycomb carbon allotropes have been 
widely investigated for this purpose following the discov- 
ery of grapheneii^; hexagonal boron nitride is another 
example^i^. With increasing chemical unit cell complex- 
ity, hexagonal transition metal dichalcogenides have been 
exfoliated into atomically thin films using mechanical 
transfer and liquid-phase sonication^^— ; many of these 
materials are semiconductors in both bulk and mono- 
layer phases, making them ideal for use in single-sheet 
nanoscale optoelectronic o^i^d"' ^ ^ . In this work we explore 
the next generation in the family of two-dimensional 
atomic crystals: stoichiometric layers of gallium chalco- 
genides (Ga2X2). 

In the bulk form, GaS and GaSe are indirect-gap 
semiconductors^, with the latter being well-known for 
its nonlinear optical properties'^, while GaTe is a direct- 
gap semiconductor used for electrothermal threshold 
switchingi^. All three materials exhibit the same struc- 
ture: they are layered compounds of Ga2X2 stoichiome- 
try in which each layer consists of two AA-stacked hexag- 
onal sublayers of gallium atoms sandwiched between two 
hexagonal sublayers of chalcogen atoms (X), as illus- 
trated in Fig. [TJ In a bulk material, these layers are 
bound in a three-dimensional structure by van der Waals 
interactions. 

In this work we discuss the optical and electronic 
properties of single-layer crystals of Ga2S2, Ga2Se2, and 
Ga2Te2. We present the electronic band structures and 
optical absorption spectra for all three materials. In 
each case we find that the band gap is indirect, with the 
valence-band maximum shifted from the T point (an in- 
verted Mexican-hat dispersion) and the conduction-band 
minimum (CBM) located at the M point in Ga2S2 and 
Ga2Te2, and at the F point in Ga2Se2. For all three mate- 





FIG. 1. (Color online) Structure of monolayer gallium chalco- 
genides Ga2X2 (X=S, Se, Te) in top and side views. 



rials we find an unusual, slightly anisotropic Mexican-hat 
shape in the electron dispersion near the edge of the va- 
lence band. However, the strongest peak in the optical 
absorption spectra — dominated by a transition between 
two bands with states which are even with respect to 
z — 5- — z symmetry — is not caused by the Van Hove sin- 
gularity associated with the Mexican-hat spectra, but is 
provided by two, near-parallel branches in the dispersion 
of the conduction bands near the Brillouin-zone corners 
K and K'. We also find that hole-doped Ga2X2 under- 
goes a Lifshitz transitio»i^ — a change in the topology of 
the Fermi surface — at hole concentrations ns = 7.96-10^'^ 
cm-^ nse = 6.13-1013 ^m-^, and nxe ^ S.Sd-lO^^ ^nr-^. 
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TABLE I. Structural parameters (as defined in Fig. [T| of 
single-layer Ga2X2 (in A) from the DFT calculations. 



X 




jLDA jLDA 
"Ga-Ga "X-X 




■PBE 
"Ga-Ga 


,PBE 

"x-x 


s 


3.53 


2.41 4.57 


3.63 


2.47 


4.64 


Se 


3.71 


2.40 4.73 


3.82 


2.47 


4.82 


Te 


4.01 


2.40 4.92 


4.14 


2.47 


5.02 



TABLE II. The band gaps of Ga2X2 according to the HSE06 
functional {Ahse) and the spin-orbit splitting in the LDA 
band structure at the K point in the highest valence (AE^'q ) 
and lowest conduction {AE^'q) bands. 



X Ahse 






S 3.28 eV 


19 meV 


21 meV 


Se 2.69 eV 


2 meV 


36 meV 


Te 2.10 eV 


14 meV 


102 meV 



II. RESULTS AND DISCUSSION 

The band-structure analysis of single-layer Ga2X2 
leading to the above conclusions was performed using 
dcnsity-functional-theory (DFT) as implemented in the 
CASTEpi^ and VASpii plane-wavc-basis codes. Tests 
have shown that the two codes yield nearly indistin- 
guishable results for the materials studied here. To 
calculate the geometries, phonon dispersions, and op- 
tical absorption spectra, we used semilocal exchange- 
correlation functionals: the local density approxima- 
tion (LDA) and the Perdew-Burke-Ernzerhofi^ (PBE) 
functionals. The screened Heyd-Scuseria-Ernzerhof 06 
(HSE06) functionaliS was used to obtain the electronic 
band structures to compensate at least partially for the 
underestimation of the band gap by semilocal function- 
als. The plane-wave cutoff energy was 600 eV. During 
relaxations a 12 x 12 Monkhorst-Pack k-point grid was 
used, while band structures were obtained with a 24 x 24 
grid. The optical absorption spectra were obtained with 
a very dense grid of 95 x 95 k-points. The artificial out- 
of-plane periodicity of the monolayer was set to 20 A. 
Phonon dispersions were calculated in CASTEP^S using 
the method of finite displacements in a 4 x 4 supercell, 
with 11 x 11 k-points, a 408 eV plane- wave cutoff, and 
an artificial out-of-plane periodicity of 15.88 A. 

Full geometry optimization was performed using both 
the LDA and PBE functionals. Table |T] shows that the 
lattice constants increase with the atomic number of the 
chalcogcn atom, while the Ga-Ga bond lengths hardly 
change. The bond lengths obtained with the PBE func- 
tional are systematically larger than those optimized with 
the LDA, as expected^!. 

The calculated electronic band structures are summa- 
rized in Fig. [2] All three materials are indirect-gap semi- 
conductors, primarily due to the valence-band maximum 
(VBM) lying somewhat off the P point. It is possible to fit 
an inverted Mexican-hat polynomial to the valence-band 
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FIG. 2. (Color online) HSE06 band structures (solid red lines) 
for Ga2S2 (top), Ga2Se2 (middle), and Ga2Te2 (bottom). The 
zero of energy is taken to be the Fermi level and the bottom 
of the conduction band is marked with a horizontal line. For 
comparison, the semilocal band structures are also shown. In 
the case of the LDA the two symmetric bands that give rise 
to the main feature in the optical absorption spectrum are 
marked in bold. 
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TABLE III. The efTective masses {m* /me) of Ga2X2 at the 
high-symmetry points in the conduction band according to 
the HSE06 functional (in units of electron mass). 



X r= K'^ 


M%r 




S 0.23 0.62 


1.52 


0.29 


Se 0.17 0.56 


3.03 


0.23 


Te 0.14 0.47 


0.43 


0.19 



dispersions Evb around the valence band maximum: 

3 

iJys = ^^2»fc'' + £^6fc®cos6^, (1) 

1=0 

where k and Lp are the polar coordinates of the wave vec- 
tors measured from the T point in units of 1/A, (p is 
measured from the F — K line, and the energy is in units 
of cV. The coefficients are listed in Table |Vl This fit 
should provide a good starting point for a simple analyt- 
ical model of the valence band in these materials. Note 
however that this fit is designed to describe the immedi- 
ate vicinity of the VBM and the saddle point, and it no 
longer describes the T point correctly, overestimating the 
r point valence band energy by 52 meV, 70 meV, and 32 
meV in Ga2S2, Ga2Se2, and Ga2Te2, respectively. 

Further analysis of the valence band reveals a saddle 
point along the F — M line, illustrated in Fig. [31 This 
saddle point gives rise to a Van Hove singularity in the 
density of states quite close to the Fermi level, which is 
found to be at the valence band edge. Due to the presence 
of these saddle points, hole-doping causes Ga2X2 to un- 
dergo a Lifshitz transition when the hole concentration 
reaches the critical value where all states are depleted 
above the energy of the saddle point, since this leads to a 
change in the topology of the Fermi surface. The carrier 
density where the Lifshitz transition takes place in each 
material is listed in the last column of Table |V] and was 
obtained by integrating the DFT density of states from 
the saddle point to the valence band edge. 

We find that the conduction-band minimum (CBM) 
is at the M point in Ga2S2 and Ga2Te2, and at the F 
point in Ga2Se2. The HSE06 band gaps of Ga2X2 are 
summarized in Table |TT1 The HSE06 band gap is ex- 
pected to underestimate the quasiparticle band gap by 
no more than 10 and is known to be applicable to 
two-dimensional materials as well2^. Our finding that 
the gap is indirect is in agreement with a recent DFT re- 
sult obtained for single-layer Ga2Se2^, although density- 
functional tight-binding calculations^^ disagree with our 
results for Ga2S2. The effective masses at the high- 
symmetry points in the conduction band are summarized 
in Table IIIII The effective mass is isotropic at F and K, 
but not at the M point. 

The semilocal band structures are also plotted in Fig. [5] 
for comparison. The LDA and PBE functionals give very 
similar results to the HSE06 functional up to the Fermi 
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FIG. 3. (Color online) Energy contours of GaX in the valence 
band using the fitted formula of Eq. ((l|. The separation of 
contours is 2 meV in all figures and the contour corresponding 
to the energy of the saddle point along the F — M line is 
highlighted. 
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level, but above that significant discrepancies arise. This 
is most notable in the case of Ga2S2, where the position 
of the CBM is ambiguous: the LDA predicts that the 
CBM is at the M point, in agreement with HSE06, while 
the PBE functional puts it at the F point. In Ga2Se2 
and Ga2Te2 both the LDA and the PBE hmctionals pre- 
dict that the CBM is at the same place as in the HSE06 
calculation. 

In the case of the semilocal DFT calculations we also 
took spin-orbit coupling into account using a relativistic 
DFT approachii. As can be seen in Fig. |3] some of the 
bands exhibit spin splitting, including the highest valence 
(AE^'q) and lowest conduction (AE^q) bands near the 
K point, most significant in Ga2Te2 (see Table lll| . While 
we were unable to calculate the spin-orbit splittings in 
HSE06 due to limited computational resources, we expect 
that they will exhibit a similar magnitude to that found 
in the semilocal band structures. 

We have calculated the phonon dispersions for Ga2X2 
using the force-constant approach in up to 4 x 4 super- 
cclls. The results are depicted in Fig. [51 suggesting that 
isolated atomic crystals of gallium chalcogenides, Ga2X2, 
are dynamically stable. 

Table IIVI shows an orbital composition analysis of the 
LDA bands of Ga2S2 around the Fermi level. We have 
found that these bands are dominated by s- and p-typc 
orbitals. Furthermore, states in each band are either odd 
or even with respect to z — > —z symmetry (this informa- 
tion is obtained from the complex phases of the orbital 
decomposition in VASP). Therefore, the interband ab- 
sorption selection rules require that photons polarized in 
the plane of the two-dimensional crystal are absorbed by 
transitions between bands with the same z — >■ — z symme- 
try of wave- functions (even— >even, odd— )-odd), and pho- 
tons polarized along the z axis cause transitions between 
bands with opposite symmetry (even— ?>odd, odd— >cvcn). 
Note that in the case of Ga2Tc2 these numbers only pro- 
vide a qualitative description of whether the orbitals are 
of mostly s 01 p character due to the strong spin-orbit 
coupling. 

The calculated optical absorption spectra (using 
semilocal DFT) arc shown in Fig. [B] The intensities 
were normalized by using graphene as a benchmark. We 
calculated the absorption of graphene using the same 
method in the 0.8-1.5 eV energy range where monolayer 
graphene absorbs 2.3 % of light, and used the calcu- 
lated graphene absorption intensity to obtain a scaling 
factor that scales the calculated intensities to the exper- 
imental value. Since in the case of Ga2S2 we found dis- 
agreement between the PBE and the IISE06 regarding 
the location of the CBM, we restrict ourselves to the 
LDA in the following. Note that the LDA results are 
only qualitatively accurate and should only be used for 
a comparative study of the different Ga2X2 monolayers 
and for an order-of-magnitude estimate of the expected 
peak positions. A better description would require a 
computationally much more expensive calculation using 
the GW approximation and the Bethe-Salpeter equation 
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FIG. 4. (Color online) Comparison between the LDA band 
structures of Ga2X2 with and without spin-orbit coupling 
(SOC). 
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TABLE IV. Orbital analysis of the LDA bands around the Fermi level in Ga2S2 at the F and K points. The labels Q, A, 
and correspond to the bands labeled as such in Fig. [21 The data were obtained by projecting the orbitals in the plane-wave 
basis set of VASP onto spherical harmonics. Dominant contributions were found to originate from s and p type orbitals; + — 
subscripts refer to even (+) and odd (— ) states with respect to z — >■ —z reflection. 



Band 


F 


K 


0+ 


0.014s*^='+0.058pz°''+0.004s^+0.189p2^ 


0.058s<3='+0.169pz°''+0.044p^PyS 


A_ 


0.142s°='+0.008pz°'-+0.098s^+0.056pz^ 


0.190s°='+0.008pz'^^+0.051p^Py^ 




0.075s°='+0.107pz*^='+0.082s^+0.001pz^ 


0.030p:,Py°"+0.036p:,p/ 



TABLE V. CoefHcients E2i (in units of eVA^') for the in- 
verted Mexican-hat dispersion near the VBM in Eq. ([l}. 
The root mean square of the residuals in Ga2S2, Ga2Se2, 
and Ga2Te2 was found to be 0.48, 0.68, and 0.67 meV re- 
spectively. The fitting parameters were obtained by fitting 
the formula of Eq. fl} to the DFT data within the range of 
0.28 < |k| < 0.42, 0.22 < |k| < 0.36, and 0.12 < |k| < 0.26 in 
Ga2S2, Ga2Se2, and Ga2Te2, respectively (with k in units of 
A~^). The zero of energy is set to the valence-band maximum. 
The last column shows the critical hole concentration {nx , in 
units of cm~^) where the Lifshitz transition takes place (see 
text). 



X 


Eo 


E2 


E4 


Es 


Eg 


nx 


s 


-0.086 


1.38 


-6.8 


5.6 


1.634 


7.96 ■ 10^^ 


Se 


-0.059 


1.49 


-11.4 


17 


4.49 


6.13 • 10^^ 


Te 


-0.044 


2.50 


-47 


240 


30.1 


3.54 • 10^^ 



III. CONCLUSION 

We have shown using first-principles density- 
functional-theory that tvifo-dimensional Ga2X2 (X=S, 
Se, and Te) crystals are stable indirect-band-gap 
semiconductors with an unusual inverted Mexican-hat 
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for excitonic correctional. In Fig. |6] we show the en- 
ergy dependence of the absorption coefBcients of various 
stoichiometric Ga2X2 monolayers. The absorption spec- 
trum starts with a shoulder at low energies, which orig- 
inates from the vicinity of the M and F points. In all 
cases the in-plane absorption is much more significant 
and is dominated by peaks at around 5 eV/4 eV/3.5 eV 
for X=:S/Sc/Tc. Analysis of the band structure reveals 
that the bulk of this spectrum comes from transitions 
between the valence band and the 2 — s- — z even band 
above the conduction band around the K point where 
the two bands are near-parallel (this is the only part in 
the band structure where the energy difference between 
a filled and an empty band matches the energy of the 
main peak). In that range, the absorption coefficients of 
Ga2X2 are comparable to and even exceed that of mono- 
layer and bilaycr graphcnc; therefore we suggest that ul- 
trathin films of GaX biased in vertical tunneling transis- 
tors with graphene electrodes could be used as an active 
element for detection of ultraviolet photons. 



FIG. 5. (Color online) LDA phonon dispersion curves for 
Ga2S2, Ga2Se2, and Ga2Te2. 



valence band. The presence of saddle points along the 
F — M line leads to a Lifshitz transition in the event of 
hole doping for which we have calculated the critical 
carrier density. We have provided an analytical fit of the 
valence-band edge and given a qualitative description 
of the optical absorption spectra, which suggest that 
ultrathin films of GaX biased in vertical tunneling 
transistors with graphene electrodes could be used as an 
active element for detection of ultraviolet photons. 
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FIG. 6. (Color online) Absorption coefBcient (the imagi- 
nary part of the dielectric function e) of various Ga2X2 two- 
dimensional crystals normalized to absolute units after it was 
compared to Im(e) evaluated for graphene in the range 0.8-1.5 
eV, where monolayer graphene absorbs 2.3 % of light (see text 
for details). Vertical lines mark the full width at half max- 
imum for the main peak. The insets illustrate the relevant 
part of k-space for the transitions between bands marked by 
circle and diamond in Fig. [2]which give rise to the main peak 
(light red circle) and the low-energy shoulder (black circle). 
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